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Abstract 

We present the results on the sohton excitations in carbon nanotubes (CNT) 
using Brenner's many-body potential. Our numerical simulations demon- 
strate high soliton stability in (10, 10) CNT. The interactions of solitons and 
solitary excitations with CNT defects are found to be inelastic if the excita- 
tions and defects length scales are comparable, resulting in a substantial part 
of soliton energy being distributed inhomogeneously over the defect bonds. 
In these solitary-cap collisions the local energy of few bonds in the cap can 
exceeds the average energy by an order of magnitude and more. This phe- 
nomenon denoted as "Tsunami effect" can contribute dynamically to the re- 
cently proposed "kinky chemistry" . We also present results of changes in the 
local density of states and variations in the atomic partial charges estimated 
at different time instants of the solitary Tsunami at the nanotube cap. 
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I. INTRODUCTION 



Since the discovery of carbon nanotubes by lijima, [|I| a rich variety of carbon nanotube 
morphologies have been experimentally observed. Carbon nanotubes (CNT) are very intrigu- 
ing objects for both experimental and theoretical investigations due to the many unusual 
properties they exhibit. Their potential for practical applications, especially in the area of 
nanotechnology, is very promising. Among their many uses envisioned include applications 
in nanoscale electronic devices, field-emission displays and quantum wires. [0 Nanotube 
doping and structure modification result in the formation of heterojunctions, diodes, 
quantum dots, P field effect transistors ^ and one-electron conductors. 

The dynamical properties of carbon nanotubes are of great interest due to their po- 
tential for useful practical applications. |§] The theoretical methods for explaining these 
properties range from accurate ah initio methods to approximate empirical schemes. Few 
of experimental findings, however, were analyzed using empirical methods that take into 
account harmonic contributions such as effects of different vibrational modes, explicitly. In 
particular, nonlinear localized excitations can transfer energy and be involved in various 
processes of interest. There have been some speculations on the role played by solitary 
excitations in heat transfer, polymer destruction and other processes occurring in molecular 
systems. While in most cases of interest harmonic approximation is sufficient to describe 
molecular processes, some molecular materials contain atoms interacting through nonlinear 
inter-atomic potentials which can give rise to soliton excitations. The investigation of these 
nonlinear effects is, therefore, interesting and timely. 

In this work we investigate the nonlinear effects in carbon nanotubes giving rise to solitons 
using empirical methods. We employ Brenner's nonlocal many-body potential for carbon 
systems. [|10[ Its simple analytical form as well as its reliability for carbon systems coupled 
with its computational efficiency makes it a natural candidate for such investigations. The 
nanotube chosen for our simulations is a (10,10) tube containing 1000-11000 atoms. The 
simulations were performed for free, rigid and capped tube ends. Our investigations include 



2 



soliton stability, reflection from free and capped CNT ends, solitons collisions as well as the 
influence of thermal fluctuations. We find that the interaction of solitary excitations with 
the CNT cap can be either elastic or inelastic. In the latter case some part of the energy is 
retained as an internal excitation of the cap with the remaining portion distributed highly 
inhomogeneously over the C-C bonds. This effect is found to be most pronounced for the 
non symmetrical caps. The degree of energy concentration depends sensitively on the cap 
structure and solitary energy and profile. In some cases the energy of few bonds in the cap 
exceeds the average solitary energy by an order of magnitude or more. Similar phenomenon 
is also observed for other defects in CNTs. 

The paper is organized as follows: In Sec. |I| we derive a nonlinear expansion of the 



Brenner's potential for longitudinal and radial atomic displacements. Sec. |TT| contains 
the formulation of an analytical equation for the longitudinal displacements of atoms along 
the CNT axis. In this section we derive the one-dimensional Korteweg - de Vries equa- 
tion assuming angular homogeneity and neglecting the radial displacements and obtain the 
corresponding supersonic soliton solution. Results of the numerical simulations of soliton 
behavior in CNTs are given in Sec. Sec. |V] contains detailed analysis of solitary colli- 
sions with CNT defects. The possible soliton contributions to various chemical and physical 
properties of CNT are briefly discussed in Sec. ^with special attention paid to the changes 
in the local density of states and variations in partial atomic charges at different instances 
of the soliton-cap interaction. Along with summarizing our results, we also propose possible 
ways of soliton generation in this section. 



II. NONLINEAR EXPANSION OF THE BRENNER'S POTENTIAL 

Carbon nanotubes (only single-walled CNT are considered in this work) are highly 
anisotropic and ordered objects with great aspect ratio. The CNT surface is formed by 
a graphitic sheet folded into a cylinder with bond lengths and angles differing slightly from 
graphite on account of the strain induced by the folding. The inter-atomic C-C potential is 
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highly unharmonic and can be approximated using classical many-body potentials contain- 
ing either exponential or power terms. Generally, these potentials do not provide an exactly 
solvable dynamical equation in the nonlinear continuum approximation and, therefore, some 
approximations needs to be made in order to yield desired solutions. 

The unharmonic solutions in nonlinear ID lattices are frequently approximated by the 



Korteweg - de Vries (KdV) equation. |T2| Since CNTs are quasi-lD objects, in order to 
obtain nonlinear excitations in CNT, it is reasonable to approximate CNT as a ID lattice 
with nonlinear effective potential. We then make the continuum approximation for this 
system and derive the corresponding KdV equation. 

We obtain the nonlinear effective potential for CNT by expanding Brenner's many body- 
potential in Taylor series up to the third order terms. The corresponding Newtonian equa- 
tions of motion in finite differences are then derived and finally, the continuous analogue 
(nonlinear partial differential equation) for the system is obtained. This results in the KdV 
equation which describes longitudinal nonlinear excitations (solitons) in CNT. 

According to Brenner's many-body potential, [|I^] the bond energy between adjacent 
atoms i and j is expressed as, 

4 = i^f-fi.,v^/. (1) 

where and are, respectively, exponential repulsive and attractive terms: V^{rij) = 
27.27 exp[-3.28{rij - 1.39)], ^^(rij) = 33.27 exp[-2.69 (r^j - 1.39)]. Bij represents an en- 
vironment dependent many-body coupling between atoms i and j containing geometric 
information associated with the system. |[T^ The total energy of the system is obtained 



by summing Eqn. |I| over all bonds. The energy and distances are measured in eV and 
Angstroms, respectively. 

The reliability of Brenner's potential has been demonstrated in many numerical calcula- 
tions of carbon systems in different phases (graphite, diamond, fuUerenes and large clusters), 
in both ground state as well as in non-equilibrium states including chemical reactions. |TI| 
Furthermore, the simple form of the potential allows us to easily derive analytical expres- 



sions for forces. This also allows efficient large scale simulations of the system using available 
computer resources. 

In the present work we consider only armchair (m, m) CNTs. The procedure and results 
are similar in the case of zigzag (m, 0) nanotubes. We use the cylindrical coordinate system 
{R, Z) for obvious reasons and align the tube axis along z-axis (Fig. |I|). With this choice 
there are 2 m atoms in the CNT layer having the same 2;-coordinates. The equilibrium 
distance between layers, io, has a slight dependence on the CNT diameter (see Table 1). 
For a (10,10) tube, io ^ 1.26A. 

Let us consider a cylindrically symmetrical disturbance in the CNT geometry in which all 
2m atoms in the n-th layer have identical z- and radial displacements from their equilibrium 
positions and -R°: Zn = Z^ + and Rn = + Pn, where Z„ and Rn are the perturbed 
coordinates and (n and p„ are the displacements from the equilibrium. Then, the coordinates 
of 2-th atom in n- th layer are (/2° + pn, , Z^ + (n). If C and p are much smaller when 
compared to characteristic length scales (Cn <^ and pn <^ -R°), then the interatomic 
potential (Eqn. p can be expanded in the Taylor series and the terms up to the third order 
retained: 

E = Eo + Y. En, (2) 

n=l 

where Eq is the ground state energy of the relaxed CNT and the perturbation energy is 
the sum over all layers of CNT. For every layer En = 2mE'^ is the sum of perturbation 
energies of 2 m atoms in the n-th layer due to atomic displacements, and energy of an atom 
E^ = 1 /2 { El + E2 + E^) , where Ei , E2 and E^ are energies of bonds emerging from any 
atom in the n-th layer (see Fig. |l|). It is convenient to reduce both ( and p to dimensionless 

units: (n Cn/4, Pn Pn/^0- 

Then an expansion E'^ over atomic ( and p displacements in the n-th layer has a form: 
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-^n — (Cn+1 — CnY — (Cn+1 ~ Cn)^ 

^1 2 ^2 / \2 ^3 / \2 / \ ^4 S 

+ ^ (Cn+1 ~ Cn) (Pn+1 + Pn) ~ ^ (Cn+1 — Cn)^ (Pn+1 + Pn) + y (Cn+1 " Cn) (Pn+1 " Pnf ■ 

(3) 

Numerical values for the expansion coefficients in Eqn. ^ are given in Table I for nan- 
otubes (5, 5), (10, 10), (15, 15) and (20,20). All coefficients in Eqn. ^ are measured in eV 
and variables C, and p are chosen to be dimensionless. Here only the leading terms in the 
nearest-neighbor interaction are retained. In this assumption the CNT is treated as a one- 
dimensional lattice with 2 m carbon atoms at each lattice site interacting through an effective 
potential given by Eqn. |^. 

Terms in Eqn. ^ can be categorized in three groups. The first contains only longitudinal 
variables the second contains only the radial variables p, and the third contains cross 
terms. The coefficients ai,a2 in the ffist group have the largest absolute values and are 
mostly independent on the CNT diameter (see Table I). The coefficients in other groups 
{hi and q) have smaller numerical values but strong dependence on the CNT diameter: the 
larger the CNT diameter, the smaller the values of these coefficients. This indicates that the 
contribution from the longitudinal degree of freedom to the dynamical properties of CNT is 
dominant, and the role of other members decreases with the increase in the CNT diameters. 
We, therefore, consider only the longitudinal degree of freedom in the study of nonlinear 
dynamics. Analytical and numerical results for the case in which both longitudinal and 
radial degrees of freedom are considered will follow in a future report. We note that these 
results differ only in minor qualitative details, the main features being determined by the 
longitudinal degree of freedom. 

Note that the oscillatory equation for the radial vibrations 

gives the lowest-frequency "breathing" mode in harmonic radial approximation for CNTs, 



and the Raman frequency for this mode for the (10, 10) CNT is ~ 165 cm~^, where M is 
the mass of carbon atom. 

The Newtonian equation of motion for longitudinal (^^displacements (neglecting radial 
degrees of freedom) is: 



M 



[ (Cn+l — Cn) — (Cn ~ Cn-l) ] " ^2 [ (Cn+1 ~ Cn)^ ~ {Cn — Cn-1 



(4) 



This is a nonlinear equation in finite differences. In the next section we derive its continuum 
analog and examine the consequences. 

III. KORTEWEG DE VRIES SOLITONS IN CNT 



Eqn. ^ presents a true ID nonlinear problem with only one variable (. The corresponding 
continuum equation is obtained using the work of M.Toda. [|l^ Let the relative displacement 
of all atoms in n-th and {n — l)-th layers be given by Xn = Cn ~ Cn-i- The continuum 
approximation is then valid if only the excitations with wave lengths much greater than the 
lattice constant Eq are allowed. The expression for the ID expansion of Brenner's potential 
(Eqn. 0) in Taylor series can be written as: 



The corresponding equation of motions for relative displacements is: 

d^Xn 



(5) 



0-2 Xn - «1 Xn 



'^2 Xn+l ~ 0-1 Xn+1 



0-2 Xn-1 ~ Xn-1 



(6) 



Eqn. ^ contains terms with Xn±i and Xn±i- is convenient to introduce a shift operator 
/ d\ 

— exp ±— — for a shift along the discrete chain by one step to the right (left) such 
\ on J 

that V^^ f{n) = /(n± 1), where f{n) is an arbitrary function of discrete coordinate n. Eqn. 
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can then be written as: 

d'^Xn 



M- 



(«2 Xl 



ai Xr 



(7) 



or 



2 sinh 



2dn 



(^2 xl - ai Xn) 



(8) 



resulting in, 



X =0. 



(9) 



If sinh ^2^^ third term in Eqn. |^ is expanded in Taylor series up to the linear 



term, then the equation can be rewritten as 



2 smh - -— 
I 2 an, 



^2 

X + 02 ^ = 0, 



(10) 



and can be factored into the following products 

a2 d 
Tiidn 



(11) 

Eqns. ( p!0| ) and (|TT]) are equivalent modulo terms of third order. Eqn. (|TT|) has a form of 



a product of two operators Fi F2 X- If "we make a further assumption: F2 x = 0, [|T3[ then 
the wave packet moving to the right (left) is described by 



M-±2V^smh 



(12) 



1 d 



If sinh ^2^J ( |T^) is further expanded in Taylor series up to the third order terms and 
n is substituted by continuum variable z, one gets the corresponding continuum equation: 

| + fa^ + fef|-%x|^=0. (13) 

at oz oz-^ oz 



where ki = k, k2 = k/24, k^ = — k, and k = ±Jai/M. 



Eqn.(p!3D can be reduced to the standard KdV equation {ur + Quu^ + u^xx = 0) by the 
linear substitution of variables: 



5 = k/24:, T = St, X = z — kt, u = a2kx/QSai. 



(14) 



Then soliton solution of flTsD has a form 



X = ±Asech'^{B z±Ct). (15) 
Here A = l2{k2/h)B'^, C = B {ki + 4^2 52). Solution of Eqn. [T| is a soliton with 



amplitude A of either compression {A < 0) or elongation {A > 0) propagating in the positive 
or negative direction along the z-axis with velocity Vgoi = C/B. It describes a continuum 
set of one-parametric solitons, i.e., coefficients A and C can be expressed through the single 
parameter, e.g., soliton half-width w ~ 1.76/ B. 

Soliton velocity can be expressed as Vsoi = A;i + 4A;2-B^ = f sound (1 + 0.52 /w"^), where 
V sound = ki = yJai/M - is the longitudinal sound velocity (~ 20 km/sec). This value is in 
good agreement with the value 20.35 km/sec obtained in Ref. fl^. The solution (Eqn. |15|) is, 



thus, the supersonic soliton. For the solitons used extensively in our numerical simulations 
(Sec. |V|), the value w = 2 was usually employed, and Vsoi ~ lASvsound- Note, that the 
solution (Eqn. |15]) has a form of solitary wave for relative ^^-displacements. For an absolute 
2;-coordinates it transforms into a kink solution ^ tanh{B z±C t). In the limit w oo the 
soliton velocity tends to the sound velocity. The dependence of Vsoi on w is shown in Fig. 
Ij. Solitons with w < 1 (less than the lattice constant io) are highly unstable in discrete 
systems and, therefore, this dependence is shown only for w > 1.5. 

The soliton in Eqn. ^ contains both kinetic and potential energies and their contribution 
to the total energy of soliton (i?*°*) are roughly the same. The dependence of (£^*°*) on w as 
well as individual contributions from kinetic {E^^^) and potential {E^°^) energies are shown 
in Fig. |[ 



The solution (Eqn. [T5|) was obtained using rather crude approximations, neglecting 
higher order terms in the expansion of Eqn. |^ and radial displacements. In order to obtain 
a more realistic soliton behavior, we next perform detailed numerical simulations. 

IV. NUMERICAL SIMULATION OF SOLITONS IN CNT 



We present our results of the molecular dynamical simulation using Brenner's many- 
body potential for the soliton evolution in the (10, 10) CNT subject to different boundary 
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conditions that include free, rigid or capped ends. The findings for the (6, 6) CNT are 
essentially the same. 

Our CNT consisted of 1000 — 11000 carbon atoms (number of layers Ni = 50 — 550), 
with a nanotube length in the range 63 — 700 A. The integration MD step was 0.35 fs and 
Verlet-Beemann algorithm was used to integrate the equations of motion. Initial conditions 
were chosen in the form given by Eqn. |l^ with typical range of parameters values: w = 
2-5, A = 0.2-0.03 A. 

Solution (Eqn. |l^) describes the compression or elongation soliton with variable ampli- 
tudes and velocities depending on the half-width w, propagating on the CNT surface. This 
solitary excitation is homogeneous with respect to the azimuthal angle $. As a first step 
in the numerical simulation we investigate the soliton evolution in CNT with frozen radial 
displacements, i.e. by allowing only the ^-displacements of the atoms. 

The time evolution of a single soliton propagating along the (10, 10) nanotube is shown 
in Fig. ^. The soliton is rather stable as it traverses the CNT length (^ 70 nm) with only 
a minimal decrease in amplitude. For CNT with cyclic boundary conditions the soliton can 
travel much longer paths without suffering any noticeable decrease in amplitude. The 
value, approximately proportional to the elastic soliton energy in the harmonic approxima- 
tion, is plotted in Fig. |[ The wider is the soliton, the greater its stability. In the limit 
of infinite CNT diameter, Eqn. |l^ describes a soliton propagating on the graphite surface. 
Its behavior is likely to be similar to the soliton in an isotropic hexagonal lattice with the 
Lennard- Jones inter-atomic potential, propagating along the [110] direction. |p!5| 

While the soliton with half-width w < 2 interacts with open CNT ends inelastically, a 
wider soliton reflects from the open ends without any appreciable changes. Rigid CNT ends, 
on the other hand, are found to preserve the soliton stability for any w. 

The collision of two solitons with different parameters is illustrated in Fig. |^. As seen in 
the figure, the solitons collide and pass through each other without changing their profiles 
and velocities, demonstrating their characteristic soliton feature. 

Our numerical simulations performed at T = 300 K showed no evidence of the infiuence 
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of thermal fluctuations on the soliton stability. Also, there were no detectable soliton distor- 
tions in the presence of defects (point mass defects, vibrationally excited bonds etc.) in the 
CNT structure. Furthermore, no significant interaction with low-intensity Raman-active 
tangential mode (~ 1585 cm~^) was observed. Results of the MD simulations, therefore, 
support the fact that the KdV equation is a good approximation for the description of 
longitudinal nonlinear excitations - solitons in CNTs - if radial displacements are frozen. 

We have also performed MD simulations of soliton evolution (propagation, reflection 
and collision between two solitons) in a ID lattice (a true ID analogue of CNT) consisting 
of atoms interacting through the potential described by Eqn. ^. Our results are in full 
agreement with the findings for the CNTs, suggesting that the neglect of higher order terms 
in the expansion of Brenner's many-body potential (Eqn. |l|) is justified. The situation 
changes, however, if the radial degrees of freedom are allowed. The folding of the graphite 
sheet to form CNT results in the mixing of longitudinal and radial modes. As a result, 
equations of motion for these two modes now contain cross-terms, which can influence the 
behavior of a longitudinal soliton. 

The numerical simulation of the initial state evolution chosen in the form of longitudinal 
soliton in CNT (Eqn. |T^) with all allowed degrees of freedom shows the soliton to be less 
stable, with a slow dissipation of the initial energy through the harmonic radial oscillations. 
Results of analytical analysis and detailed numerical investigations of the longitudinal-radial 
solitons will be presented in a forthcoming paper. 

V. INTERACTION OF SOLITARY EXCITATIONS WITH NANOTUBE CAPS. 

We next investigate the interaction of the solitons with nanotube caps. Our results for 
the interaction of narrow {w = 2 — A) solitons with non-local defects (free tube ends and 
caps) shows it to be inelastic with the excitation energy being distributed approximately 
evenly over all the bonds in this region. Interesting effects can be expected in the case when 
the length scales of both the excitations and defects are comparable. Towards this end, we 
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perform a detailed numerical investigation of interaction of excitations with the CNT caps. 

Solitons with w = 10 carry too small an energy to induce any noticeable effects in the 
cap. So we used solitary excitations with other profiles. When a solitary excitation has a 
characteristic length scale much greater than the lattice constant the nonlinear effects can be 
ignored and can be treated within the harmonic approximation. Furthermore, for long-wave 
excitations, the dispersion effects in harmonic approximation are not significant. In this 
case, an arbitrary excitation propagates along the harmonic chain with speeds close to the 
sound velocity with only a minimal change in its profile. We have examined the propagation 
of a longitudinal solitary excitation using the parameters x = Asech^[B [z — Vsoundi)] with 
w — 10 and A — 0.1 A. This excitation was found to be highly stable even when all degrees 
of freedom were allowed and survived the reflection from free ends. 

Same results were obtained for solitary waves with other analytical profiles and similar 
widths and amplitudes. The other types of excitations studied included an excitation of the 

sound ^) 



Gaussian profile x = ~j= ^^V 



, also with w — 10 and ^4 = 0.1 A. 



2^2 

Although the soliton propagation is accompanied by a highly coherent and predominantly 
longitudinal displacements of the atoms on the cylindrical surface of the CNT, the atoms on 
the cap behave quite differently after the sohtary wave - cap collision. When the solitary 
excitation reaches the cap, there is an accumulation of energy at the cap. Part of the 
total energy of the solitary wave is likely to be dissipated and/or accumulated on structural 
defects. We consider an example in which solitary wave interacts with non symmetrical 
caps. 

Our simulation of soliton interaction with the nanotube cap consisted of the following 
sequences: (A) The incoming excitation causes the cap to "infiate" with the cap attaining 
nearly a spherical shape. At this instant the potential energy of the cap attains a maximum 
value, although the energy distribution in the bonds are unequal due to the asymmetry in 
the cap. Atoms in the tip of the cap are displaced by ~ 0.3 A in the positive z-direction; 
(B) The recoil of the elastic energy of the cap causes atoms to move in the opposite direction 
and at some instant the kinetic energy attains a maximum value. The velocities of atoms 
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are all predominantly directed along the z-axis although they differ significantly from each 
other in their magnitudes; (C) Atoms at the cap continue to be propelled by inertia and the 
elastic potential energy attains a maximum again. At this instant, the cap acquires regions 
with negative curvature, and the displacement of the tip of the cap from its equilibrium 
position is ~ —0.8 A; (D) After this the atoms in the cap start to move in the opposite 
direction toward the equilibrium state and second, but less intense, maximum in the kinetic 
energy is observed. 

The nanotube cap continues to oscillate with rapidly decreasing amplitude, the energy 
being dissipated in heating the cap and the CNT "body" . Note the large total amplitude of 
the tip displacements (> 1.0 A). A striking feature of this preliminary results observed in 
MD simulations is the inhomogeneity of the energy distribution in the cap after the solitary 
wave - cap collision. We next investigate the cap dynamics and energy distribution over 
atoms in more details. 

A typical time dependence of the kinetic energy distribution over atoms in the (10, 10) 
cap is shown in Fig. ^ as a histogram. Schlegel diagram of this cap is shown in Fig. 

This cap was chosen from a full list of 9 342 topologically different caps with isolated 
pentagons generated previously. [|T^] Atom numbering is common to both Figs. ^ and ^b, 
and atoms with numbers less than 30 (in the cap base) are not shown as they have lower 
(< 0.02 eV) excitation levels. Two "waves" of excitations are visible in this Figure: one at 
t = 750 — 800 MD steps (stage (B)) and the other with lower excitation level - at 1300 — 1600 
MD steps (stage (D)). The atoms with energy > 0.1 eV are drawn in large circles in Fig. |^a; 
atoms with energy > 0.05 eV - in medium circles. An excitation area is localized and bears 
a larger portion of total cap energy. One can see that the kinetic energy is concentrated 
mainly at the tip of the cap and irrespective of the location of pentagons. 

We denote this effect of energy concentration as "Tsunami effect" (named for the effect 
explaining ocean waves coming to a beach) because initially long-wave and low-amplitude 
excitation is concentrated into a sharp impulse of energy if the conditions for the solitary 
wave propagation are changed in a special manner. It was found that this effect depends on 

13 



the cap structure and is more pronounced in less symmetrical caps. The effect is observed if 
the cap radius is comparable with solitary wave half-width w (in the case shown in Fig. |^ for 
w = 10). Otherwise, the Tsunami effect is less prominent. Analogously similar phenomena 
with varying degrees of perfections were observed for other defects (kinks and bends) in 
CNT. 

We have also calculated electronic properties of the capped CNTs during the solitary 
wave - cap collision using a generalized tight binding scheme of Menon and Subbaswamy. 



17| The variation of the partial Mulliken charge distribution obtained using this scheme is 



shown in Fig. |^. These calculations were performed for time instants corresponding to the 
maximal cap distortions (stages (A) and (C)). The most significant change, as expected, is 
observed for atoms at the tip of the cap (indicated by large number lables in Fig. |^) where 
perturbation in the structure is a maximum. 

The local electron density of states (DOS) averaged over atoms at the cap at stages (A) 
through (C) of solitary wave - cap collision are shown in Fig. |^. The DOS is obtained 
from the Green's function G(E) = [ES-H+i5]^^, where H is the generalized tight binding 
Hamiltonian, S the overlap |T^, and 6 = 0.05 eV. Note the rather large variations in the 
DOS. 



VI. CONCLUSION 

We have, thus, demonstrated analytically and using numerical simulations that the KdV 
solitons are reasonably good approximations for the description of nonlinear excitations in 
CNTs. The soliton stability increases with the CNT diameter and, in the limit R oo, the 
solution (Eqn. |T5|) describes the nonlinear excitations in graphite. Our numerical simulations 
of the evolution of the initial longitudinal soliton confirms its high stability. 

We have also performed numerical simulations of the interaction of solitary excitations 
(different from solitons) with defects in CNTs, especially with caps (the Tsunami effect). Our 
simulations show many new fascinating features in the dynamics of CNTs. The collision is 
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found to be partially inelastic, and the cap excitation is highly inhomogeneous. Interestingly, 
at some instances the energy in a few bonds considerably exceeds the averaged energy of 
solitary excitation. This process can provide extra dynamical contribution to the recently 
discovered phenomenon of the enhanced reactivity of defect sites in CNTs termed "kinky 
chemistry". [|TB],|Tg 

The issues concerning ways of solitons excitation in real CNTs and their possible con- 
tribution to various chemical and physical phenomena are not yet clear. In reality, the 
generation of the highly coherent ultra short longitudinal displacements of atoms in the 
CNT seems hardly probable. This, however, is not the case for the long wave solitary exci- 
tations. In this connection we would like to point out recent results on the opto-mechanical 
effect in CNT. [pO| , pl| Conceivably a short flash of light could trigger the CNT excitation in 



such a manner that the solitary excitations will be generated. They can also be generated 
due to external factors such as electron or ion impacts, stress release and other mechanisms. 
If solitary excitations do exist in CNTs in reasonable "concentration", then they can be 
identified through their contribution to some detectable CNT properties. A possible exam- 
ple is the Tsunami effect when the energy, smoothly distributed in the large scale excitation 
of an arbitrary profile, can be concentrated on a few bonds of a defect with considerable 
energy excess. The greater the distortions in CNTs, the greater is this influence. It can also 
be a way for a "self-healing" of CNTs when a rearrangement of structural imperfection is 
activated by solitary excitations. This effect also can promote chemical reactions such as C2 
incorporation into nanotube caps. 



The high specific heat of a rope of SWNT observed [22| can be partially explained by 



solitary excitations generated in parallel with phonons. Heat transfer in CNTs is another 



example. p3|-p5|] Emission of short electric field pulses, when solitary wave inelastically 
interacts with the non-symmetrical caps or other defect sites, can be yet another example. 
Our results of the calculation of variation in charge distribution and DOS confirm large 
deviation of the electronic properties from their equilibrium values at the solitary excitation 
- cap collision. The process of charge distribution variation under dynamical excitation is 
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probably inverse to the opto-mechanical phenomenon in CNTs. pO[, Nanotechnology is 
one area where the mechano-electric property of CNT can be used. 

We had initially planned to title our paper "Solitons in carbon nanotubes" , but after the 
preparation of this manuscript we became aware of a very recent work on soliton by Chamon 
2^ with the same title. In this paper the spontaneous lattice distortions are investigated 



similar to the case of polyacetylene. The solitons are topological domain walls separating 
different symmetry-broken vacua with different Kekule bond-alternation structures. Note 
that the solitons discussed by Chamon and those found in the present work have different 
nature: we have considered solitons formed due to the elastic and non-linear properties of 
CNT, while Chamon solution is a topological soliton. 
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FIGURES 




Z„.|Z„ Z„,| 



FIG. 1. A part of (10, 10) nanotube. The tube is oriented along the z-axis. Atoms in 
n, (n lb l)-th layers are marked by black and gray circles, respectively. 



Soliton half-width (w) 

FIG. 2. The dependence of soliton velocity Vsoi on the soliton half-width w. Vgound is the 



longitudinal sound velocity. 
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FIG. 3. The dependence of total soliton energy on soliton half-width w and individual 



contributions from kinetic and potential (E^"^) energies. 








FIG. 4. Soliton evolution in the (10, 10) nanotube with 3 000 atoms (150 layers) with free ends. 
The soliton half-width isw = 2 and = Xn is the square of relative longitudinal displacements of 
the n-th and (n — l)-th layers. The modulation of the soliton amplitude in this figure is explained 
by the lattice discreteness and soliton narrowness. Actual x values are shown in insertion at two 
different time instants (labeled by full circle and empty stars, respectively). Only maximal values 
of |Ap reflect the soliton stability. 
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FIG. 5. Temporary evolution of two solitons with w\ = 2 (right at t = 0) and ii)2 = 3 (left at 
t = 0) moving initially toward each other. CNT and notations are the same as in Fig. 
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(b) 




FIG. 6. (a) Schlegel diagram and the atomic numbering of non-symmetrical (10,10) CNT cap; 
(b) temporary evolution of the kinetic energy distribution over atoms in the cap. Initial condition 
- solitary excitation with w = W moves from the tube to the cap and at t ~ 600 MD step reaches 
the cap. The maximum energy concentration is attained at t » 800 MD steps. 
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FIG. 7. Variation of the partial Mulliken charges at the cap at different time instants vs. 
charges of the relaxed cap. Atomic numbering is the same as in Fig.^ (A) - largest "positive" 
expansion of the cap (stage (A)); (C) - "negative" shrinkage (stage (C)). 
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FIG. 8. Variation of the densities of states averaged over the cap atoms at different instances. 
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TABLES 



TABLE L Coefficients (in eV) in the expansion of Brenner's potential (Eqn. ^ and the 
equiUbrium distance between layers io (in A) for CNTs of different diameters. 







fin in"! 


{10, 10) 




ai 


1 no on 
iUz.oU 


102. /6 


1 no 7n 
102. ^^0 


1 no ^it; 
102. DO 








^aO.OV 




Ol 


1 1 sn 

il.OU 


9 QS 


1 "^9 


u. / U 




U.o4 


0.10 


n n /I 
0.04 


n nn 
0.00 




9.06 


4.59 


3.09 


2.31 


64 


18.42 


2.40 


0.72 


0.30 


Cl 


12.24 


6.16 


4.14 


3.10 


C2 


71.34 


35.82 


24.06 


17.94 


C3 


152.40 


154.56 


154.20 


153.66 


io 


1.2601 


1.2573 


1.2568 


1.2565 



24 



